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ABSTRACT 

We present an XMM-Newton detection of two low radio surface brightness SNRs, G85.4+0.7 and 
G85.9— 0.6, discovered with the Canadian Galactic Plane Survey (CGPS). High-resolution XMM- 
Newton images revealing the morphology of the diffuse emission, as well as discrete point sources, are 
presented and correlated with radio and Chandra images. The new data also permit a spectroscopic 
analysis of the diffuse emission regions, and a spectroscopic and timing analysis of the point sources. 
Distances have been determined from H I and CO data to be 3.5 ± 1.0 kpc for SNR G85.4+0.7 and 
4.8 ± 1.6 kpc for SNR G85.9-0.6. 

The SNR G85.4+0.7 is found to have a temperature of - 12 - 13 MK and a 0.5-2.5 keV luminosity 
of ~ 1 — 4 x 10 33 Z?§ 5 erg/s (where -D3.5 is the distance in units of 3.5 kpc), with an electron density 
n e of ~ 0.07 — 0.16(/I?3.5) _1 / 2 cm~ 3 (where / is the volume filling factor), and a shock age of 
- 9 - 49(/L> 3 .5) 1/2 kyr. 

The SNR G85.9-0.6 is found to have a temperature of - 15 - 19 MK and a 0.5-2.5 keV luminosity 
of ~ 1 — 4 x 10 34 £> 2 g erg/s (where -D4.8 is the distance in units of 4.8 kpc), with an electron density 
n e of - 0.04 - 0.10(/£> 4 .8)" 1/2 cm" 3 and a shock age of - 12 - 42(/L> 4 .8) 1/2 kyr. 

Based on the data presented here, none of the point sources appears to be the neutron star associated 
with either SNR. 

Subject headings: ISM: supernova remnant — ISM: individual object: G85.4+0.7 — ISM: individual 
object: G85.9-0.6 — stars: neutron — X-rays: ISM 



1. INTRODUCTION 

In 2001, two new supernova remnants (SNRs) with 
low radio surface brightness, G85.4+0.7 and G85.9— 0.6, 
were di scovered in Canad ian Galactic Plane Survey 
CCGPSldTavlor et all l2003h data and conf irmed in X- 
rays with ROSAT data (jKothes et al.|[200ll ). Both show 
distinct shells in the radio band with an extended re- 
gion of X-ray emission in the centre. The radio surface 
brightness of G85.4+0.7 at 1 GHz is £i G Hz < 1 X 10~ 22 
Watt m -2 Hz -1 sr -1 and the radio data also indicate 
that the SNR has a non-thermal shell with angular di- 
ameter « 0.4° which is surrounded by a thermal shell 
with an angular diameter of w 0.6° and is located within 
an H I bubble. The bubble also contains two B stars 
which may have been part of the same association as the 
SNR's progenitor star. G85.9— 0.6 has a radio surface 
brightness Sighz < 2 x 10" 22 Watt m~ 2 Hz -1 sr" 1 and 
it has no discernible H I features, indicating that it is 
expanding into a low density medium, perhaps between 
the local and Perseus spiral arms. The most likely event 
which would produce an SNR in such a region would be 
a Type la supernova. 

1 Department of Physics and Astronomy, University of Mani- 
toba, Winnipeg, MB, R3T 2N2, Canada. 

2 Canada Research Chair 

3 Department of Physics, George Washington University, Wash- 
ington, DC 20052, U.S.A. 

4 National Research Council of Canada, Herzberg Institute of 
Astrophysics, Dominion Radio Astrophysical Observatory, P.O. 
Box 248, Penticton, BC V2A 6J9, Canada. 

J Department of Physics and Astronomy, University of Calgary, 
2500 University Drive NW, Calgary, AB T2N 1N4, Canada. 

6 Department of Physics & Astronomy, Brandon University, 270- 
18th Street, Brandon, MB, R7A 6A9 Canada. 



X-ray observations are important to the study of SNRs, 
particularly those with low surface brightness, because 
they provide information about the morphology and 
emission processes of these objects, which are indica- 
tors of both the nature of the supernova explosion which 
formed them and the properties of the progenitor star. 
SNRs with low surface brightness are expected to be 
formed after the core collapse of a massive star in a 
Type Ib/c or Type II explosion, because the stellar wind 
would have blown away much of the interstellar medium 
(ISM) surrounding it, leaving a low ambient density into 
which the shock from the supernova expands. A Type 
la supernova can also result in a low surface brightness 
SNR if the surrounding ISM has a low density, such as it 
would if it were located between two spiral arms of the 
galaxy. Thermal X-ray emission from a SNR arises as the 
blast wave of the explosion travels through and shocks 
the ISM, and as a reverse shock travels back into and 
shocks the ejecta. The X-ray spectrum of the SNR gives 
information about the temperature, the density, and the 
luminosity of the shocked material, while imaging data 
provides information about the size and morphology of 
the region. 

The low angular and spectral resolutions as well as the 
small number of counts in the ROSAT data did not allow 
spectroscopy nor detailed imaging to be done, so XMM- 
Newton observations, described in |}2j have been made 
in order to confirm the detection of the SNRs, and to 
perform imaging, spectroscopic and timing studies. De- 
tailed X-ray imaging, described in |}3l is used to map the 
diffuse emission and compare it to the location and size 
of the radio shells, and Chandra data have been used to 
search for compact objects not resolved by XMM. Spec- 
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tral parameters obtained in §4] lead to an estimate of 
such quantities as temperature, density, and luminosity 
of the SNRs. In §5] the point sources are catalogued and 
an attempt at identification is made by matching their 
positions to objects in other catalogues. Timing analysis 
is performed to identify any pulsar candidates. The dis- 
tances to the SNRs are derived from H I and CO data in 
§6] The results are discussed in §7] 

Preliminary results w ere presen ted in 

iJackson. Safi-Harb. fc Kothesl (pool and iSafi-Harbl 
(12006ft . 

2. OBSERVATIONS 
2.1. XMM-Newton Observations 

SNR G85.4+0.7 was observed with XMM-Newton on 
2005 May 31 for 11.5 ks (Obs ID: 0307130101, PI: S. 
Safi-Harb) and again on October 27 for 15.2 ks (Obs ID 
0307130301, PI: S. Safi-Harb), because of proton flares in 
the first observation which made a large fraction of the 
data unusable. A 29 ks XMM-Newton observation was 
made of SNR G85.9-0.6 on 2005 November 24 (Obs ID 
0307130201 PI: S. Safi-Harb). 

The PN (jStriider et al.ll2001h and MOS ijTurner et all 
2001) data were reduced with the latest version of SAS 
(6.5.0) and events during proton flares were filtered out 
for producing images and spectra by using the SAS rou- 
tines evselect to create rate files and tabgtigen to gener- 
ate good time intervals (GTIs) for filtering the event files 
with the evselect routine. This rendered the first obser- 
vation of G85.4+0.7 useful only for imaging, and slightly 
reduced the integration time, from 29 ks to 26 ks for 
the G85.9— 0.6 observation. The second observation of 
G85.4+0.7 shows no evidence for proton flaring and the 
removal of proton flares did not significantly reduce the 
integration time. Images and spectra were created from 
the event files using evselect. The total effective exposure 
time for G85.4+0.7 was 14 ks for PN and 16 ks for the 
MOS instruments and for G85.9-0.6 it was 23 ks for PN 
and 26 ks for MOS. 

To facilitate source detection, exposure maps were 
made using eexpmap and the images were combined us- 
ing emosaic. The spectra were binned using a minimum 
of 25-50 counts per bin for the MOS instruments and 50 
or 100 counts per bin for the PN for the point source 
and diffuse spectra respectively, using grppha. This lat- 
ter grouping was necessary because the background sub- 
traction added excessive noise to the spectra, particu- 
larly at high energies, where the spectra are background- 
dominated. The background for both SNRs was calcu- 
lated using various methods described in §4.11 

2.2. Chandra Observations 

A 14.5 ks observation of SNR G85.4+0.7 was made 
on 2003 January 26 (Obs ID: 3898, PI: S. Safi-Harb) by 
the Advanced CCD Imaging Spectrometer (ACIS-S; G. 
Garmire 7 ). SNR G85.9-0.6 was observed with ACIS-S 
for 20.2 ks on 2003 October 3 (Obs ID: 3899, PI: S. Safi- 
Harb). Both observations were made at a focal plane 
temperature of — 120°C. 

The analysis of the Chandra data was done using 
Chandra Interactive Analysis of Observations (CIAO 



version 3.3 s ). For both Chandra observations, the data 
were corrected for charge transfer inefficiency (CTI) 
with tools provided b y the ACIS team of Pennsylvania 
(|Townslev et al.ll2000D . Events with ASCA grades (0, 2, 
3, 4, 6) were retained, and periods of high background 
rates were removed. Data from hot pixels were elimi- 
nated. The effective exposure time for the observation of 
SNR G85.4+0.7 was 14.3 ks and for SNR G85.9-0.6 it 
was 19.9 ks. The backgrounds for the point source spec- 
tra were extracted from annuli surrounding the sources. 
The spectra were grouped with a mimimum of 30 counts 
in each bin. Data from the S2 and S3 chips were used. 
Data from the other chips were not used because they do 
not overlap the XMM field of view. 

3. IMAGING 

The previous X-ray images from the ROSAT All-Sky 
Survey did not resolve the SNRs and the point sources in 
the field, so the XMM-Newton PN and MOS data can be 
used to determine the nature of the sources. The results 
are described below. 

3.1. G854 +0.7 

The 0.5-2.5 keV XMM-Newton PN and MOS mosaic 
image of G85.4+0.7 with radio contours overlaid is shown 
in Figure [TJ The image includes only energies between 
0.5 and 2.5 keV in order to match the ROSAT bandpass 
and because there is negligible diffuse emission above 2.5 
keV. The X-ray image is smoothed with a Gaussian fil- 
ter with a 3 pixel radius. The X-ray emission is in the 
same location and is a similar angular siz e as indicated 
by the ROSAT data, shown in Figure 6 of iKothes et alJ 
(|2001h . with an approximate angular radius of 6.7'. It is 
much better resolved, however, and the regions of diffuse 
emission and the point sources are clearly visible. 

The top two panels of Figure [2] show the 0.5-2.5 keV 
and 2.5-10 keV XMM-Newton images of G85.4-0.7, 
smoothed in the same way as Figure [TJ The absence of 
diffuse emission in the hard X-ray band is clear from the 
middle panel. The lower panel shows the 0.5-2.5 keV 
image with sources subtracted, smoothed further and 
scaled to emphasize the diffuse emission, and smoothed 
contours are overlaid to show the morphology of the rem- 
nant. The contours indicate a centrally filled morphology 
with a slightly elongated shape. The circles in the top 
panel indicate spectral extraction regions described in 
gHand gH 

There are nine point sources in the image from which 
spectra can be extracted, and the analysis is described 
in SJSJ The extraction regions for the spectra, including 
the background for the diffuse emission, are indicated 
and labeled in Figure [2j and the spectral analysis of the 
diffuse emission is described in §321 

The Chandra image of SNR G85.4+0.7 is shown in 
Figure [3l The extraction regions used for the source 
and background regions are shown, and the applicable 
regions from the XMM-Newton observation are overlaid. 
The wavdetect tool in CIAO is used for source detection. 
Six of the point sources identified with XMM-Newton (1, 
2, 3, 4, 6 and 7) are clearly visible on the Chandra im- 
age. Twelve additional sources appear on the Chandra 
image (labeled 10—21) but the spectra extracted from 
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these point sources do not contain enough counts to en- 
able meaningful spectral analysis, so while they are po- 
tential neutron star candidates and sources 10, 11, 12, 
13, 15, and 16 are situated well within the radio shell 
of the SNR, further investigation of these sources is not 
possible at this time. 

3.2. G85.9-0.6 

The XMM-Newton PN and MOS mosaic image of 
G85.9— 0.6 with radio contours overlaid is shown in Fig- 
ure [H The image is smoothed in the same way as the 
G85.4+0.7 image, with a Gaussian filter with a three- 
pixel radius. As with G85.4+0.7, the angular radius 
of approximately 5.8' and location of the X-ray emis- 
sion shown in this figure matches that of Figure 6 of 
iKothes et al.l (|2001l ) and many features are resolved. 0.5- 
2.5 keV and 2.5-10.0 keV X-ray images are shown in the 
top two panels of Figure [5l It is clear from the hard 
X-ray image that the SNR does not appear above 2.5 
keV, and the X-ray spectrum is background dominated 
above 2.5 keV, as will be explored in §4.31 The bottom 
panel of Figure O shows the same data as in the top panel 
with point source emission subtracted, scaled to empha- 
size the diffuse emission. The contours in the bottom 
panel are smoothed to show the overall morphology of 
the remnant. The diffuse X-ray emitting region seems 
to be round in shape, indicating a centrally filled X-ray 
morphology. The circles in the top panel indicate spec- 
tral extraction regions described in £14.31 and £17.41 

Spectra were extracted from 8 point sources and the 
diffuse region, as indicated in the top panel of Figure 
The analysis of the spectrum from the diffuse region is 
described in £|4.3i and S}5] describes the analysis of the 
point sources. 

Figure [6] shows the Chandra image of SNR G85.9— 0.6. 
The locations of the extraction circles for sources 1 and 
8 in the XMM-Newton image are shown, in addition to 
the additional point sources found on the Chandra image 
(labeled 9—27). The Chandra image does not show any 
point sources which have sufficient counts to do indepen- 
dent spectroscopic analysis, and only one of the sources 
identified from the XMM-Newton image (source 8) lies 
within the field of view of the Chandra image, with which 
a Chandra source has been identified. 

4. SPECTROSCOPY 

The spectra from the regions of diffuse emission were 
extracted from areas indicated in Figures[2]and[5]for both 
PN and MOS. The diffuse emission regions and point 
sources are indicated and labeled in the Figures. 

The spectra of the diffuse emission for both SNRs 
were fit to various models: a simple bresstrahlung 
model with gaussian lines, VMEKAL, VNEI, and VP- 
SHOCK, each modified by interstellar absorption, in 
this cas e the Wisconsin absorpt i on mo del (wabs in 
XSPEC, iMorrison fc McCammonl (fl983T) ). VMEKAL 
is a collisional ionizat ion equilibrium mo del based on the 
model calculations of iMewe et al.l (|1985h of the emission 
spectrum from hot diffuse gas with Fe calculations by 
iLiedahl et alj (|1990t ). with variable abundances. VP- 
SHOCK is a non-equilibrium ionization (NEI) model 
with variable abundances which models a plane paral- 
lel shock in a plasma with constant electron temperature 
T and a range of ionization timescales with an upper 



limit r = n e t, where n e is the postshock electr on density 
and t is the shock age (|Borkowski et al.|[200ll ). VNEI is 
a NEI model similar to VPSHOCK, except with a single 
ionization timescale. 

4.1. Background Estimation 

Spectra of diffuse faint regions require particular at- 
tention to be paid to the choice of background region. 
Instrumental emission lines can dominate the spectrum if 
the background region is poorly chosen, leading to errors 
in spectral fits and a discrepancy between the PN and 
MOS spectra. The intuitive background region would be 
a relatively source free one at approximately the same 
galactic latitude as the source to minimize contamina- 
tion by the Galactic ridge. However, because the effect 
of instrumental emission lines is not uniform with po- 
sition on the CCDs, particularly on the PN instrument 
(|Lumbl I2002T ) , and exceeds the variation with galactic 
latitude for these observations, the background regions 
were chosen as closely as possible to the point opposite 
the source region through the centre of the image, tak- 
ing care that there is no overlap between the source and 
background regions and omitting regions around point 
sources. The background regions for both SNRs are nec- 
essarily smaller than the source regions, but very little 
change in the spectral fits resulted from choosing vari- 
ous background regions in the same general area of the 
image. 

An attempt has been made to estimate the MOS back- 
ground spectra using the XMM-Newto n Extended Source 
Analy sis Software (XMM-ESAS) (jSnowden fc Kuntzl 
2006). As yet, background estimation for the PN instru- 
ment is not included in this package. Data from filter 
wheel closed observations and from the unexposed cor- 
ners of the MOS CCDs are used to generate background 
spectra. This serves as a test of the backgrounds gener- 
ated from the observation itself, described above. To fit 
spectra using the XMM-ESAS backgrounds, it is neces- 
sary to include s ome additional spe c tral c omponents, as 
recommended in lSnowden fc Kun tz (2006). Unabsorbed 
Gaussian lines of zero width at energies of 1.49 and 1.75 
keV are fit to the MOS spectra in addition to the model 
describing the emission from the SNR. The 1.75 keV line 
was also added to the fit to the MOS spectra using the 
observation background. It was not necessary to add 
the 1.49 keV line to the fit of the G85.4+0.7 MOS spec- 
tra using observation backgrounds because it was ad- 
equately subtracted, but the 1.49 keV line was added 
to the G85.9— 0.6 MOS spectra using observation back- 
grounds. It was also necessary to include a low energy 
broken power law component to the G85.9— 0.6 fit of the 
ESAS-subtracted spectrum, as instructed in the ESAS 
documentation, to correct for instrumental effects, be- 
cause the statistics allowed residuals at low energies to be 
clearly seen in the ESAS background subtracted spectra. 
The broken power law was not added to the spectrum of 
G85.4+0.7 because the spectra were noisy and did not 
require an additional model component. Results of spec- 
tral fits using both background estimation methods for 
each SNR are given in Table [T] The parameters in the 
left hand column of Table Q] for each of the SNRs are 
from the simultaneous fit to the MOS spectra using the 
ESAS background subtraction method and the PN spec- 
trum with the background spectrum extracted from the 
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observation itself, and the right hand column contains 
the parameters from the simultaneous fit to the MOS 
and PN spectra from both of which has been subtracted 
the observation background. 

As a third background subtraction method, we have 
attempted to use blank sky event files that are avail- 
able for the PN and MOS instruments and generated 
to enable the production of background spectra for ex- 
tended sources. These event files comprise a superpo- 
sition of pointed observations from which sources have 
been removed. Background files produced from blank sky 
event files therefore simulate the detector response in an 
actual observation, including any instrumental emission 
lines. Events are first selected from the original event 
files based on sky p osition using the SelectR ADec script. 
The sky cast script (|Read fc Ponmanll2003f) is then used 
to cast the new event file onto the sky coordinates for 
the particular observation. The background spectra can 
be extracted from these files from the same region as 
for the source spectra. The BACKSCAL keyword is ad- 
justed in the background file to scale it to the source 
file. Unfortunately, the resulting background subtracted 
spectra are oversubtracted for the PN instrument, and 
contain negative counts below 0.6 keV and above 1.5 
keV. Thus the blank sky background spectra cannot be 
used for this analysis. Possibly the oversubtraction re- 
sults from the fact that there are no observations within 
30° of the SNRs G85.4+0.7 or G85.9-0.6 contained in 
the blank sky datasets, and therefore systematic errors 
arise from the different level of background emission. The 
background subtracted PN spectra of G85.4+0.6 using 
a background from the observation itself and from the 
blank sky file is shown in Figure El The spectra from 
which the blank sky background was subtracted were not 
used and thus the parameters do not appear in Table [TJ 

4.2. G854 +0.7 

The PN and MOS spectra between 0.5 and 2.5 keV are 
simultaneously fit in order to determine the best param- 
eters while reducing instrument-specific and systematic 
effects. The upper limit of 2.5 keV is used for the spec- 
tral fits because the spectra are background dominated 
above this energy as shown in Figure [TJ 

Fitting to the VMEKAL model results in a reduced x 2 
value of 1.76 for the ESAS background or 1.63 for the 
observation background (this refers to the MOS spectra 
only; the background spectra extracted from the obser- 
vation was used for every PN spectrum), and a kT value 
of 0.65 ± 0.03 or 0.65 ±0.01 keV (all errors are 2a) for 
the ESAS and observation background respectively. 

The spectra are best fit with an absorbed VPSHOCK 
with with kT of 1.1 t° i keV for the ESAS back- 
ground or 1.0^0 2 keV for the backgrounds extracted 
from the observation itself, and ionization timescale of 
(6.4±|'|) x 10 10 cm" 3 s for the ESAS background or 

(8.0^34) x 10 10 cm" 3 s for the background extracted 
from the observation. The spectrum for the diffuse re- 
gion is shown with the ESAS background in Figure[8]and 
with the observation background in Figure [5J The pa- 
rameters for the simultaneous PN and MOS fits, with 2a 
uncertainties, are given in Table [TJ along with derived 
quantities such as luminosity and age. The fitted pa- 
rameters agree well within uncertainty for the two back- 



ground estimation methods. The abundances are for the 
most part consistent with solar, with exceptions given 
in Table [TJ An analysis of the spectral fit parameters 
has shown G85.4+0.7 to have a 0.5-2.5 keV luminos- 
ity of (3.1±i j) x 10 33 erg s" 1 or (2.0±^) x 10 33 erg s" 1 
for the ESAS or observation backgrounds respectively, in 
both cases based on the estimated distance of 3.5 ± 1.0 
kpc for G85.4+0.7, determined in Sj6j The normaliza- 
tions of the instrumental lines in the MOS spectra from 
which the ESAS background has been subtracted are 
(1.7±0.2)xl0~ 5 and (5.5±?;g) x 10~ 5 photons cm^s- 1 
for the lines at 1.49 and 1.75 keV respectively. The nor- 
malization for the line at 1.75 keV for the observation 
background is (3.7+}'5) x 10~ 5 photons cm _2 s _1 . 

To verify the fits to elemental abundances, the VNEI 
model was used in place of VPSHOCK. All elemental 
abundances were found to be consistent within error to 
be solar or subsolar, except O and Fe. O was clearly 
enhanced above a solar value, giving strength to the VP- 
SHOCK result. The elemental abundances using the 
VMEKAL model were also qualitatively similar to the 
VPSHOCK result. 

4.3. G85.9-0.6 

The approach used for fitting the diffuse spectra of 
G85.9-0.6 is similar to that used for G85.4+0.7. Again, 
the upper limit of 2.5 keV is used for the spectral fits be- 
cause the spectra are background dominated above this 
energy. The PN diffuse spectrum was fitted simultane- 
ously with the MOS spectra from which either the ESAS 
or observation background was subtracted. 

Fitting to the VMEKAL model results in a reduced 
X 2 of 2.21 for the ESAS background and 2.48 for the 
observation background, a kT value of 0.70 ± 0.01 or 
0.68 ± 0.01 keV. 

The spectra of the diffuse emission from the remnants 
are best fit with an absorbed VPSHOCK model with 
kT of 1.3±^ keV for the ESAS background or 1.6±g;| 
keV for the background extracted from the observation 
and ionization timescale (6.8^q g) x 10 10 cm -3 s for the 
ESAS background or (5.1±^) x 10 10 cm" 3 s for the 
background extracted from the observation. The fit- 
ted parameters with 2a uncertainties are given in Ta- 
ble [TJ and the diffuse spectrum with the ESAS back- 
ground is shown in Figure [10] and with the observa- 
tion background in Figure 1111 The abundances are for 
the most part consistent with solar, and the exceptions 
are shown in Table [TJ The luminosity of G85.9— 0.6 
is (2.4±*; 2 ) x 10 34 erg s^ 1 for the ESAS background 

or (2.7^2 3) x 10 34 er S s_1 f° r the observation back- 
ground, in both cases based on the estimated distance 
of 4.8 ± 1.6 kpc for G85.9-0.6, determined in $3 The 
normalizations for the 1.49 and 1.75 keV lines in the 
ESAS subtracted MOS spectra are (1.22±0.08) x 10~ 4 
and (3. 7lQg) x 10 -5 photons cm _2 s _1 respectively. The 
broken power law which was added to the fit of the 
MOS ESAS-background subtracted spectrum had Ti = 
1.2t2;$, T 2 = 3.3t°; 4 , E hleak = 0.91 ± 0.05 keV, and 
a normalization of (8.1+*' 4 ) x 10~ 4 photons/keVcm 2 - 
s at 1 keV. The normalization of the 1.75 keV line for 
the observation background is (3.2 ± 1.1) xlO -5 photons 
cm s . 

When the VNEI model was used in place of VPSHOCK 
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to check the elemental abundances, as was done with 
G85.4+0.7, all elemental abundances were found to be 
consistent within error to be solar or subsolar, except O 
and Fe, both of which are clearly above solar abundance 
when the error bar is taken into account, in agreement 
with the VPSHOCK result. Using the VMEKAL model, 
the elemental abundances are again qualitatively similar 
to those obtained from the VPSHOCK model, except Mg 
is above solar. 

5. POINT SOURCE ANALYSIS 

The previous ROSAT data of the SNR regions pro- 
duced contours of the diffuse emission but did not re- 
solve any point sources. The presently considered XMM- 
Newton and Chandra data allow for point sources to be 
resolved and located within the field of view. Nine clearly 
distinguishable point sources are seen on the image of 
G85.4+0.7 (Figure©, of which six (sources 1, 2, 3, 4, 6 
and 7) are in common with point sources in the Chan- 
dra observation, as shown in Figure[3] Twelve additional 
point sources are clearly resolved in the Chandra obser- 
vation (labeled as sources 10 through 21 in Figure©, 
but they do not possess sufficient counts to allow for 
meaningful spectral analy sis to be don e . It s hould be 
noted that the B stars in iKothes et al.l (|2001h are out- 
side of the field of view of the X-ray observations. Eight 
point sources are seen on the XMM-Newton image of 
G85.9— 0.6 (Figure© of which only one source (source 8) 
is visible on the Chandra image, which is pointed toward 
the region encompassing sources 1 and 8 in Figure 
Figure [HI also shows 19 additional point sources from the 
Chandra observation (labeled 9 through 27), but none of 
them possess sufficient counts for their spectral parame- 
ters to be sufficiently constrained. In the present work, 
the nine point sources detected with XMM-Newton in 
the field of G85.4+0.7 and eight in the field of G85.9-0.6 
are analysed using spectral and timing techniques, to lo- 
cate any candidates for identification as neutron stars 
or pulsars which would have formed at the time of the 
supernova explosion. 

5.1. Source Identification 

To identify the X-ray point sources in the G85.4+0.7 
and G85.9— 0.6 fields, ewavelet, a wavelet detection algo- 
rithm which is part of the SAS 6.5 package, is used. For 
each source, the output of the routine gives position on 
the image and in sky coordinates, source counts, source 
extent, as well as errors in those quantities. Because the 
PN and MOS instruments have slightly different fields 
of view, only those sources found in the combined image 
which are also found on the PN image are analysed. The 
sources must also have an extent (size on the image) sim- 
ilar to the PSF at that location on the image, and they 
must contain enough counts (> 50) for the spectral and 
timing analysis. 

A catalogue of newly discovered X-ray point sources 
is given in Tables [2] and H The letters XMMU in the 
object designations indicate that they were discovered 
in XMM-Newton data and a prefix of CXO indicates a 
discovery in Chandra data. The ROSAT All Sky Sur- 
vey catalogue was checked to see if any of the sources 
appear there, and source 1 in the G85.9— 0.6 image (Fig- 
ure © is well within the 24" la error circle of the co- 
ordinates of the ROSAT All Sky Survey object 1RXS 



J20591 1.3+444730, and source 1 in the G85.4+0.7 image 
(Figure© lies marginally within the 24" la error circle 
of the RASS object 1RXS J205058.7+452135. No other 
point sources in this study match any in the RASS cat- 
alogue, indicating that these point sources are too faint 
to be included in the RASS catalogue. 

The sky coordinates of sources meeting the above 
criteria are searched within the extent, which is ap- 
proximately the PSF size, in vari ous catalogues, in - 
cluding the USNO-A2.0 catal ogue dMonet et al.l ll998f). 
the USNO-B1.0 cat alogue dMonet et al.l l2003h. the 
SKY2000 catalogu e (iMvers et al.ll2002f l. the catalogue 
given in iGuarinod (l!992f) . and the 2MASS catalogue 
dSkrutskie etlTE boe) , to determine whether the sources 
have already been identified in another waveband. Some 
of these catalogues give blue magnitude and if this 
is converted into flux with the relation mi — m,2 — 
—2.5log(Fi/ F2) and compared with the X-ray flux, this 
ratio is one factor that can be checked to favour identifi- 
cation as a neutron star. In the above equation, mi and 
TO2 are the blue magnitudes of the object and a standard 
star and Fi and F2 are their fluxes. Typical neutron 
stars have a flux ratio in optic al to X-ray bands of ~ 10 
(|Lvne fe Graha m-Smith] I"l998l ). whereas X-ray emitting 
O- o r B-type stars such as 77 Carinae ([Corcoran et al.l 
I2000D or r Scorpii (jMewe et all l2003h have a ratio of 
~ 10 7 , which enables an identification to be made of 
neutron star candidates. The optical to X-ray flux ratios 
for the point sources, matched with each of the catalogue 
objects within their extents are shown in Tables [2] and 
[31 These values represent the optical to X-ray flux ratio 
should the catalogue object match the X-ray source. If 
the optical catalogue source and the X-ray point source 
are not the same object, the ratio is meaningless. How- 
ever, an X-ray source without any optical counterparts 
within the PSF, such as source 6 or 7 from the G85.9— 0.6 
data, may be a good neutron star candidate. In other 
words, this test does not exclude X-ray point sources 
as neutron star candidates but rather indicates that a 
source emits in X-rays, and is much fainter in the optical 
waveband, and is unlikely to be a stellar object. The op- 
tical to X-ray ratios for all point sources for which the ra- 
tio can be calculated (ie all sources except 6 and 7 in the 
G85.9— 0.6 field) is > 10 11 , which is much greater than 
the ratio expected for neutron stars, indicating that these 
objects are not neutron stars, assuming that the optical 
sources are the true counterparts. Other potentially in- 
teresting objects (e.g. neutron star or AGN candidates) 
are Chandra sources 14, 16, 17, 18, 19, 21, 24, 25, and 27 
in the G85.9— 0.6 field, three of which have no counter- 
part at all (21, 24, and 27), and the rest of which have 
either only a 2MASS counterpart or a 2MASS counter- 
part which is a much better match for the position than 
the optical counterpart. These objects are shown in Fig- 
ure [6] and listed in Table [3] 

5.2. Spectral Analysis 

In addition to the optical to X-ray flux comparison de- 
scribed in §5.1| the X-ray spectra of the point sources can 
be examined to determine the likelihood that any of them 
are neutron stars. The X-ray spectrum of a rotation pow- 
ered pulsar is typically hard with a po wer law photon i n- 
dex (r) of approximately 0.5-2.1 (e.g. iGotthelJ (|2006h ). 
though there may be an additional blackbody component 
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from the neutron star surface which could dominate the 
spectrum, depending on the age and the star's magnetic 
field. Anomalous X-ray pulsars t ypically have photon 
indice s between 2.4 and 4.6 (e.g. IWoods fc Thompson! 
(2006)). The background-subtracted spectra of the point 
sources, where the background region consists of an an- 
nulus surrounding each point source, are fit to an ab- 
sorbed power law and the resulting fitted parameters are 
shown in Table H In addition, the 0.5-2.0 keV and 2.0- 
10.0 keV counts are given, which is particularly useful 
for estimating an X-ray hardness ratio when the fit to 
an absorbed power law yielded a large x 2 value, indi- 
cating an unsuitable model. Sources 1, 4, and 9 in the 
G85.4+0.7 field and sources 4, 6 and 7 in the field of 
G85.9— 0.6 have been identified as neutron star candi- 
dates. Sources 1 and 9 of G85.4+0.7 both have photon 
indices of ~ 2.5, which does not rule them out as neutron 
star candidates but they both have a high optical/X-ray 
flux ratio. Both source 4 of G85.4+0.7 and source 4 of 
G85.9— 0.6 have a relatively soft X-ray spectrum com- 
pared with typical pulsars (and the \ 2 value for the fit 
of source 4 on G85.9— 0.6 indicates that the power law 
model does not match the data well) and both have a 
high optical/X-ray flux ratio. However, it is not certain 
that the optical sources associated with these objects are 
the true counterparts. Neutron star candidates 6 and 7 in 
the G85.9— 0.6 field have a photon index typical to neu- 
tron stars and no optical counterpart. However, source 6 
has a column density (Ah) which is greater than that of 
the SNR itself, which indicates, along with its position 
relative to the SNR, that it is unlikely to be the asso- 
ciated neutron star. Source 7 has a Nn value which is 
similar to that of the SNR, but its position indicates that 
it is probably not associated with the SNR. When con- 
sidered separately, XMM-Newton and Chandra spectra 
of G85.4+0.7 point sources 1 and 3 yield spectral param- 
eters which agree with each other well within the range of 
uncertainty. The Chandra spectra of sources 2, 4, 6, and 
7 do not contain enough counts for independent spec- 
tral analysis but the Chandra spectra for these sources 
arc included in the analysis leading to entries in the ap- 
propriate rows of Table [4] and the spectral parameters 
of the combined XMM-Newton and Chandra spectra of 
these sources lie within the uncertainty of those for the 
XMM-Newton spectra alone. As an example, the XMM- 
Newton PN and MOS and Chandra spectra of G85.4+0.7 
source 1 is shown along with a fit to an absorbed power 
law in Figure fT2l 

5.3. Timing Analysis 

A timing analysis is done on the PN data for the neu- 
tron star candidates identified above to search for pul- 
sations. The timing resolution in PN full window mode 
(68.7 ms) for both G85.4+0.7 and G85. 9-0.6 allows for 
a search up to ~ 7 Hz, which would fail to identify fast 
rotation-powered pulsars, but would identify slowly ro- 
tating anomalous X-ray pulsars. The PN events file is 
first barycentre corrected. The photon arrival times for 
a region surrounding each source are used first in a fast 
Fourier transform (FFT) search to identify possible fre- 
quencies that should be investigated further. Around 
each frequency identified by the F FT search, Rayleigh 
(|Leahv. Eisner, fc~W cisskopfl 1983[) (also known as Zf), 
Z% and Z\ ( Buccheri et al.l 19831 ). and epoch folding 



searches are done. The statistical significances of any 
identified frequencies are calculated using the \ 2 prob- 
ability along with the number of degrees of freedom, 
which is 2n for the Z% search and ribins — 1 for the epoch 
folding search, where ribins is the number of bins in the 
lightcurve. In this case the number of bins used is 12 and 
20 for the two epoch folding tests performed on the data, 
chosen because those numbers of bins produce a rela- 
tively detailed lightcurve while maintaining a reasonable 
number of counts per bin. With regard to the Z^ test, 
the Z\ 1 Z\ 1 and Z\ searches are performed and would 
identify typical pulsar X-ray lightcurves, which usually 
exhibit either a broad variation or two or three peaks 
per cycle. 

The timing resolutions of the MOS instruments (which, 
in the data acquisition mode used for these observations, 
is 2.6 seconds) and Chandra (3.24 seconds) are not suf- 
ficient for meaningful timing analysis of this type to be 
done. 

Of the six neutron star candidates found from the point 
source spectra of the two observations, none were found 
to exhibit a periodic signal. However, a future dedicated 
timing search may uncover one of these sources as a pul- 
sar. 

5.4. Results of point source analysis 

Based on a combination of the optical to X-ray ratios 
(assuming the optical sources are the counterparts of the 
X-ray sources), distances from the SNR centres, spectral 
parameters, and comparisons of A^h value to that of the 
diffuse emission, none of the point sources is likely to 
be the neutron star candidate associated with the SNR. 
However, further observations of the Chandra objects 
may reveal one of them to be a neutron star or AGN. 

6. DISTANCES 

It is interesting to note that the absorbing H I 
column density derived from the X-ray spectra is 
smaller for G85.9-0.6 than it is for G85.4+0.7, even 
though G85.9— 0.6 w as predicted to be further away by 
iKothes et~al"1 (|2001h . The knowledge about the fore- 
ground H I column density gives us another option to 
constrain the distance to these objects. Both SNRs are 
believed to be located outside the solar circle. In the 
outer Galaxy the radial velocity is decreasing monotoni- 
cally with distance, independent of the Galactic rotation 
model used. Hence, we can integrate foreground atomic 
and molecular hydrogen as a function of distance, by 
integrating the spectroscopic data down to the appro- 
priate radial velocity. Neutral hydrogen d ata are from 
the C anadian Galactic Plane Survey (see iTavlor et al.l 
120031 for further information), and the 12 CO(1-0) molec- 
ular line data ar e from the Columbia CO survey of 
iDame et al.l (|1987l ). 

Since atomic hydrogen is usually optically thick we 
could not simply integrate the H I emission, as emission 
does not represent all of the hydrogen actually present 
along a Galactic line of sight. To correct for this we have 
to determine the optical depth of the H I along the line 
of sight. One way of doing this is looking at background 
point sources, preferably of extragalactic origin, which 
are very bright radio continuum emitters at 1420 MHz. If 
these sources are bright enough we can see their emission 
being absorbed by the foreground H I and use them to 
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probe the ISM along their line of sight through the whole 
Galaxy. The absorbing neutral hydrogen column density 
integrated over the velocity interval dv is then defined 
by N Rl {v)[cm' 2 } = 1.8224 x 10 18 T s (v) t{v) dv. r is 
the optical depth, which can be determined from the ab- 
sorption profile by t(v) = -ln((T on (v)—T B (v)/T bg ) + l), 
here T bg is the brightness temperature of the absorbed 
background source, T on is the continuum subtracted H I 
brightness temperature at the position of that source, 
and Tb is the brightness temperature of the absorbing 
H I cloud. Tb is represented by the average off-source 
H I brightness temperature, which is determined in an 
1' wide elliptical annulus just outside the absorbed back- 
ground source. T s is the spin temperature of the absorb- 



ing cloud, which is defined by T s (v) = Tb(v)/(1- 



-r(v) 



Equ ations for N-^ i(v ) , t{v) , and T s (v) can be found in 
e.g. iRohlfs fc Wilsonl (pool . 

To derive the complete foreground H I column we have 
to add twice the molecular hydrogen column density 
Nh 2 - This is derived over the velocity interval dv us- 
ing its relation to the CO brightness temperature Tg°, 
determined by iDame et all (|200lD : N H2 (v)[cm,- 2 } = 
1.8xl0 20 Tg°{v)dv. 

For the SNR G85.4+0.7 we used the absorption profiles 
of two background sources to determine the foreground 
atomic hydrogen column density. One source is located 
at the very centre of the remnant (Figure fl3l source 1) 
at £ = 85.31° and b = +0.66° and the other just to the 
south (Figure[H source 2) at I = 85.31° and b = +0.32°. 
For G85. 9-0.6 we used the bright source at I = 85.80° 
and b = —0.65°. To calculate the amount of foreground 
molecular hydrogen we averaged the Dame et al. (2001) 
CO data, which has a pixel size of 30' over the four pixels 
closest to the centre of the X-ray emission. The final 
combined H I column density profiles are displayed in 
Figure [131 

If we now compare the absorbing H I column density, 
which we derived from the X-ray spectra (see Tabled]), 
with the iVjj j - velocity diagrams in Figure 1131 we can 
get an estimate for the radial velocities of the SNRs. 
For the SNR G85.4+0.7 we derive a radial velocity of 
about —9 km s _1 averaging over source 1 and 2. This 
nicely confirms the radial velocity of —12 km s _1 , which 
was determined for the stellar wind bubble surrounding 
G85.4+0.7 by Kothes et al. (2001). For G85.9-0.6 no 
previous estimate of the radial velocity exists. The com- 
parison of the absorbing H I column (Table [4| with the 
A/tj j - velocity diagram in Figure [13] results in a radial 
velocity estim ate of —32 ± 6 km s " 1 . 

Previously, iKothes et alj (|2001h found a distance to 
G85.4+0.7 of about 3.8±0.6 kpc, based on the radial ve- 
locity of its host stellar wind bubble using a flat rotation 
model for the Galaxy with a Galacto-centric radius of 
8.5 kpc for the Sun at a velocity of 220 km s^ 1 around 
the Galactic centre. A distance of ~5 kpc was pre- 
dicted for G85.9-0.6. The radial velocities of G85.4+0.7 
(-12±3 km s" 1 ) & G85.9-0.6 (-32±6 km s" 1 ) indi- 
cate they are beyond the Solar circle (Rq =7.6 kpc, 
lEisenhauer et a]|l2005h . but from T b (£, v) diagrams m 
this direction they are not residents of the Perseus Spi- 
ral arm (which shows as a large H I feature extended in 
longitude, very nearly centred on —40 km s _1 ). 

A new kinematic-based distance method has been de- 



veloped by l^stiEi&^^Wnniffii ((2006). The approach 
is based on a model of the Galactic H I density distri- 
bution and velocity field, that is fitted to observations, 
rather than relying on a purely circular rotation model 
assigned to the object (as in standard kinematics). The 
model's density component is that of a warped thick disk 
of H I, with axisymmetric features, and a two-arm den- 
sity wave pattern in the disk. The velocity field compo- 
nent models a non-linear response of the gas to the den- 
sity wave (see iRoberts 1 4969; W ielen Ill979t ). producing 
shocks at the leading edge of major arms, and streaming 
motions within the disk, depending on the location of an 
object in the spiral phase pattern. This distance method 
has been shown to accurately reproduce spectrophoto- 
metric distances to H II regions throughout the second 
quadrant of the Galactic plane. 

The best fit synthetic H I profile (with velocity field 
as in Figure [14]) towards these objects shows that 
G85.4+0.7 is 3.5±1.0 kpc distant. This distance indi- 
cates G85.4+0.7 is within the Local spiral arm, a fea- 
ture in between the major Sagittarius and Perseus spiral 
arms of the Milky Way. While the velocity of G85.9-0.6 
is less certain compared to G85.4+0.7 (dete rmined by 
association with H I; see IKothes et al.l 1200 1 ; , it should 
be noted that in G85.4+0.7's case, comparison of the 
X-ray absorbing column to the total integrated hydro- 
gen nuclei (in H I and 12 CO data) gives a good veloc- 
ity estimate. Hence, we can be reasonably sure that 
—32 km s _1 is near to the true velocity of G85.9— 0.6 
as well. This places it within the Perseus arm spiral 
shock, which is 4.8±1.6 kpc in this direction (the arm's 
potential minimum is 5.0 kpc). For any velocity in the 
range — 40< vlsr < — 18 km s -1 , the distance range 
shown in the fitted velocity field (Figure fl"4]) is small, 
about 4.1-4.8 kpc. 

Although SNRs stemming from Type II events are pri- 
marily found within the arms (a s massive prog enitors are 
born mainly in the arm 's shock. IWiele"n~lll979f) . the loca- 
tion of G85.9— 0.6's within the shock does not necessarily 
vitiate its identification as a Type la. For example, Ty- 
cho's SNR is known to be Type la, but its velocity and 
distance clearly associate it with the Perseus Spiral arm. 
It is possible that the binary precursor of G85.9— 0.6 mi- 
grated into the arm along its Galactic orbit before ex- 
ploding as a Type la event. 

7. DISCUSSION 

7.1. Determination of shock age and mass of emitting 

gas 

To estimate t, the shock age, for the SNRs, the rela- 
tion t = r/n c is used, where r is the upper limit of the 
ionization timescale from the VPSHOCK model. The 
electron density n c is determined from the distance in 
cm -D cm , angular size in radians a, given by the diameter 
of the extraction region, and the relation njj = n c /1.2, 
where njj is the volume density of hydrogen within. 

Given that Norm = iqiI^""^ ) ; where Norm is the nor- 
malization of the VPSHOCK model, which is propor- 
tional to the emission measure, and the electron and 
hydrogen densities n c and tih are assumed to be uni- 
form, a volume filling factor / is employed so that n c = 
^ 2.88xio^(Noim) p rom these calculations, it is deter- 
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mined that for G85.4+0.7, n e w (0.11 ±0.03)(/ D3.5)- 1 / 2 
for the ESAS background (for the MOS instrument) or 
(0.1lto o4)(/D 3 .5)- 1 / 2 for the observation background, 
and t = (18tg 7 )(/ J D 3 .5) 1/2 kyr for the ESAS background 
or (23^24) (/D3.5) 1 / 2 kyr for the observation background. 
For G85.9-0.6, n e « (0.07t°;^)(/L» 4 .8)^ 1/2 for the 
ESAS background or (0.07±0.03)(/L>4.s)~ 1/2 for the ob- 
servation background and t = (30^ 1 3)(/Z?4.8) 1 ^ 2 kyr for 

the ESAS background or (22t^)(/D 4 . 8 ) 1/2 kyr for the 
observation background. These are larger than the age 
estimates from the radio data. 

The mass of the emitting gas is calculated based on 
spherical emitting regions with the size given by the ex- 
traction radius (given in Table [TJ , composed of 92% hy- 
drogen and 8% helium, and the relation rijj = n e /1.2 
is used as for the above calculation of t. The mass 
of the emitting gas in G85.4+0.7 is (2.8 ± l^D^M® 

for the ESAS background and (2.8^)0^ M G for 
the observation background and that for G85.9— 0.6 is 

(2.7t 1 ^)L>4 / 8 2 M for the ESAS background and (2.7 ± 
1.8)1)4 for the observation background. 

7.2. Background subtraction 

The background subtraction was problematic, partic- 
ularly for the spectrum of G85.4+0.7, which is fainter 
and required a different model to be fit to it depend- 
ing on which background region on the observation was 
used, leading to a very careful selection of the background 
extraction region, the process of which is described in 
Sj4Tj Instrumental lines, which are fortunately differ- 
ent for PN and MOS, nevertheless increased the value 
of x 2 for the combined fits, even when the background 
was chosen very carefully, and needed to be explicitly fit 
when the ESAS background was used, as described in 
£ 14.11 In addition, the soft spectra of the SNRs exhibit 
large amounts of noise in the high energy end of the spec- 
tra, allowing fits to be made only up to 2.5 keV. Future 
longer observations of these SNRs will help to resolve 
these difficulties, allow for better fits to the abundances, 
eliminate ambiguities in the spectral results, and per- 
haps allow for conclusive identifications of neutron star 
and pulsar candidates in addition to other point sources. 

7.3. G854 +0.7 

A comparison between the morphologies of the X-ray 
diffuse emission and the radio emission of G85.4+0.7 
can be seen in Figure [15] in which the X-ray point 
sources have been removed and the X-ray image has been 
smoothed to 1' to match the radio contours. The diffuse 
X-ray emission exhibits some structure, and it lies in the 
approximate centre of the radio shell. A spectrum was 
extracted from the central blob, the position and size 
of which is shown in the top panel of Figure [21 and it 
could not be adequately fit with a non-thermal power 
law model, indicating that a pulsar wind nebula origin 
for the emission can be ruled out. 

The fact that the SNR does not exhibit any limb 
brightening and the X-ray emission is mostly concen- 
trated in the centre, as projected in the plane of the sky 
and three dimcnsionally, can be interpreted as evidence 
that the X-ray emission is produced by the ejecta. The 



fact that there appears to be no emission from the swept 
up material can be explained if the remnant is evolution- 
arily young. 

With the ejecta interpretation, if it is assumed that 
the free electrons are evenly distributed in the extrac- 
tion area, the resulting ages are 18 and 23 kyr for the 
two backgrounds. In Figure [TBI it can be seen that the 
outer radio shell has an approximately constant radius, 
whereas the inner shell has a smaller radius in the ver- 
tical centre, which expands as the latitude changes, and 
this indicates that the SNR is moving either toward us 
or away from us. The two shells appear to meet near the 
bottom of the image, so the radius of the SNR can be 
approximated as that of the outer shell, which is ~ 0.3° 
on the image. The average velocity of the expanding 
SNR is then 910 or 730 km/s for the two ages. How- 
ever, the morphology of the X-ray emission is not very 
smooth, indicating that the filling factor / << 1. As- 
suming / = 0.2, the electron density n e would be 0.25 
cm -3 for either the ESAS and observation background, 
and the ejecta masses would be 1.3 M Q for either back- 
ground. The shock age would then be 8000 or 10300 
years, and the average expansion velocity would be 2150 
or 1640 km/s. For a 10 51 erg supernova explosion, the 
ejecta mass would be 22 or 37 M©, or 2.2 or 3.7 M© 
for a 10 50 erg explosion, indicating that the lower energy 
explosion would result in an ejecta mass consistent with 
a young freely expanding SNR. 

The radio continuum emission indicates that there is 
some swept up material. Assuming that the density in- 
side the stellar wind bubble was ~ 0.01 cm~ 3 before the 
explosion, the SNR would have swept up ~ 6.5 M© of 
material, which is on the order of the ejecta mass, and 
means that the SNR is in the transition between free 
expansion and Sedov expansion. 

The slightly above-solar abundance of O in the dif- 
fuse spectra reinforces the hypothesis that the super- 
nova most likely resulted from a core collapse, though 
the large error bars weaken the argument. The abun- 
dances of all elements other than O and Fe are below 
solar, but they would be enhanced if the X-rays were 
from the ejecta. However, it is possible that the spectral 
parameters for the abundances are affected by the poor 
quality of the spectra and by the uncertainties associated 
with the background subtraction. 

Since the most likely origin of G85.4+0.7 was a core 
collapse supernova (see S|6|), it is possible that a neutron 
star or pulsar which is associated with this SNR can be 
found. Given its power law X-ray spectrum, proximity 
to the centre of the remnant, and similar JVh value to 
the diffuse emission, source 1 in Figure [2] is possibly the 
associated neutron star, but sources 4 and 9 are within 
the radio shell and are also candidates, given their spec- 
tral parameters. However, the low fit quality of source 
4 (x 2 ~ 2) indicates that it is not likely to be a neu- 
tron star, and source 9 is situated far from the centre of 
the diffuse emission so is less likely to be the associated 
neutron star. Chandra sources 10, 11, 12, 13, 15, or 16 
in Figure [3] can also be considered as neutron star candi- 
dates (provided the optical counterparts listed in Table [2] 
are not the true counterparts. 

7.4. G85.9-0.6 
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Figure [TBI shows X-ray and radio images of G85.9— 0.6, 
produced in a similar way to Figure 1151 The diffuse 
X-ray emission appears to contain less structure than 
G85.4+0.7, and, as for G85.4+0.7, a spectrum extracted 
from the central blob, the position and size of which are 
shown in the top panel of Figure[5j indicates that a pulsar 
wind nebula origin for the emission can be ruled out. 

Because Figure [15] shows no limb brightening and the 
emission is mostly from the central blob, a similar argu- 
ment to that used for G85.4+0.7 can be employed here to 
suggest an ejecta interpretation for the X-ray emission. 
Assuming that the free electrons are evenly distributed 
in the extraction area the resulting age is 30 or 22 kyr 
for the two backgrounds. Using the average distance be- 
tween the centre of the X-ray emission and the shell, 
~ 0.2°, the radius of the shell is 15.3 pc The average 
expansion velocity would then be 500 or 680 km/s. The 
emission is concentrated in the inner 6', which indicates 
an electron density of 0.20 cm -3 and an ejecta mass of 
1.0 M© which agrees with the predicted mass for a type 
la explosion, 1.4 Mq. This would result in an age of 
10600 or 7800 years and average expansion velocity of 
1350 or 1850 km/s. 

For a Type la supernova, the explosion energy is ~ 10 51 
erg and the ejecta mass is 1.4 Mq. Unlike for G85.4+0.7, 
the density in the interarm region is closer to 0.1 cm , 
resulting in a swept up mass of nearly 50 Mq. This 
indicates that the SNR is in the Sedov expansion phase, 
which is described by R = l^Eo/noY^t 2 / 5 , where R 
is the radius in pc, Eq is the explosion energy in units 
of 10 51 erg, no is the ambient density in cm~ 3 , and t is 
the age in units of 10 4 years. For an explosion energy 
of 10 51 erg, a radius of 15.3 pc, and an age of 10600 
or 7800 years, the resulting ambient density is 0.84 or 
0.55 cm -3 which are consistent with the density of the 
interarm region, but the swept up mass would be 360 
or 240 Mq , from which it should be possible to measure 
thermal X-ray emission, from the part of the shell that 
is included in the X-ray pointing. The current expansion 
velocity would be dR/dt = 0A(R/t), which is 570 or 760 
km/s. 

The interpretation of G85.9— 0.6 as having been pro- 
duced by a Type la supernova is reinforced by the 
Fe abundance, which is well above solar. As with 
G85.4+0.7, the ejecta interpretation is called into ques- 
tion by the remaining abundances, which should be 
above solar, but are instead below solar. This could again 
be a result of poor quality spectra. 

Giv en that the radio results described in iKothes et al.1 
(2001), as well as the distance presented here, indicated 
that G85.9— 0.6 was most likely produced by a Type la 
supernova, it was not expected to find a neutron star as- 
sociated with this SNR. The above solar Fe abundance 
for this SNR is consistent with a Type la explosion. An 
identification of sources 6 and 7 in Figure [5] has not 
yet been made. They are bright X-ray emitting ob- 
jects with no known optical or radio counterpart, making 
them good neutron star or radio-quiet AGN candidates, 
though if one of them were a neutron star, it would be 
unlikely that it is associated with the G85.9— 0.6 SNR be- 
cause they are both situated far outside the radio shell 
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of the remnant, and furthermore, the value of Nn for 
source 6 does not match that of the SNR itself, and the 
fit quality of source 7 (\ 2 ~ 2) indicates that an ab- 
sorbed power law is not a good fit. Their identification 
with possible 2MASS counterparts (see Table [3|) makes 
them more likely radio-quiet AGN than neutron stars, 
provided the 2MASS objects are the true counterparts. 
A future detailed deep X-ray or multiwavelength study 
of these objects should be undertaken to identify and 
further study them, even though they are probably not 
associated with the G85.9— 0.6 SNR. Source 4 is on the 
edge of the radio shell of the SNR, but its identification 
as a neutron star is questionable because of its photon 
index and hardness ratio, as well as the fact that it is 
not expected that there is a neutron star associated with 
this remnant. 

7.5. Mixed Morphology Interpretation 

The centrally filled morphology of both SNRs and the 
thermal nature of their X-ray emission confined within 
the radio shells suggest that they belong to the class of 
mixed-morphology SNRs (also known as thermal com- 
posites; Rho & Petre 1997). The origin of the thermal 
X-ray emission interior to the radio shells in these SNRs 
has been attributed to several mechanisms which include: 
a) cloudlet evaporation in the SNR interior (White & 
Long 1991), b) thermal conduction smoothing out the 
temperature gradient across the SNR and enhancing the 
central density (Cox et al. 1999), c) a radiatively cooled 
rim with a hot interior (e.g. Harrus et al. 1997), and d) 
possible interaction with a nearby cloud (e.g. Safi-Harb 
et al. 2005). While modeling these SNRs in the light 
of the above mentioned models is beyond the scope of 
this paper and has to await better quality data, we can 
rule out the cloudlet evaporation model based on the low 
ambient densities inferred from our spectral fits (see §7.1 
and Table 1). Except for Fe and possibly O, the abun- 
dances inferred from our spectral fits are for the most ele- 
ments consistent with or below solar values, as observed 
in most mixed-morphology SNRs. However, enhanced 
metal abundances have been observed in younger SNRs, 
e.g. 3C 397, estimated to be ~5.3 kyr-old and proposed 
to be evolving into the mixed-morphology phase (Safi- 
Harb et al. 2005). The ages inferred for G85.4+0.7 (~8- 
10 kyr; see Table 1 and jS| and G85. 9-0.6 (-8-11 kyrs; 
see Table 1 and §7.4|) suggest a later evolutionary phase 
where only shock-heated ejecta from Fe (for G85.9— 0.6) 
and possibly Oxygen are still observed. 



XMM-Newton is an ESA science mission with instru- 
ments and contributions directly funded by ESA Mem- 
ber States and NASA. This research is supported by 
the Natural Sciences and Engineering Research Coun- 
cil of Canada (NSERC), and partly by NASA grant 
NNG05GL15G. This research has made use of NASA's 
Astrophysics Data System, and the Canadian Galac- 
tic Plane Survey, a Canadian project with international 
partners supported by NSERC. We thank an anonymous 
referee for useful comments. 
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TABLE 1 

Fits of G85.4+0.7 and G85.9-0.6 XMM-Newton PN and MOS diffuse region spectra to an absorbed VPSHOCK and 

DERIVED QUANTITIES. ERRORS ARE 2a UNCERTAINTIES. ABUNDANCES OF He, C, AND N ARE FROZEN TO SOLAR. THE Nl ABUNDANCE 

IS TIED TO THE Fe ABUNDANCE IN THE FITS. 



Parameter G85. 4+0.7 G85. 9-0.6 



MOS Background used: 


XMM-ESAS 


Region in Figure [2l 


XMM-ESAS 


Region in Figure [5] 


Nn (10 22 cm- 2 ) 


n »fi+ 04 


0.87 ±0.06 


n fis+ 03 

U - bS -0.04 


0.70 ±0.03 


kT (keV) 


i 1+0.8 
L - L -0.3 


i.o±8:3 


i n+0.5 
1 -°-0.2 


1 6 +0 - 5 


t (10 10 cm _3 s) 






b ' S -0.9 


5 1+!-5 
°' X -1.0 


Norm a 


(2.4±g:») x 10~ 3 


(2.3±°-g) x 10- 3 


(l.lt°;J) x 10- 3 


(1.1 ±0.2) x 10~ 3 


O b 


1.2 ±0.5 


1 2+ 2 ' 7 


i K+0-5 
1 -°-0.4 


i 7 +0.6 


Ne 


4+ ' 1 


5+ - 2 
u -°-0.4 


o.o±g;S 


+0 ' 2 


Mg 


0.3 ±0.1 


4+ - 1 


0.6l°; 2 


6+ 01 
U.D_ 2 


Si 


0.1 ± 0.1 


0.2 ±0.2 


1.1 ±0.2 


0.8 ±0.2 


S 


4+ - 2 




o.8±8:l 


4+ - 8 
u -*-0.4 


Fe 


i 1+0.3 
L - L -0.2 


i-i±8:l 


2 6+ ' 5 
z - D -0.2 


2.6 ±0.3 


xl (y) 


1.41 (209) 


1.38 (211) 


1.50 (460) 


1.78 (465) 


Distance (kpc) 
Radius (pc) c 


3.5 ± 1.0 


4.8 ±1.6 


(6.4±1.8)£>3. 5 


(7.5 ± 


2.5)D 4 . 8 


Temperature (MK) 


13+1 


12±I 




19±i 


0.5-2.5 keV 










absorbed fiux d 


1.1 x 10" 12 


9.6 x 10" 13 


1.5 x 10" 12 


1.6 x 10" 12 


0.5-2.5 keV 










unabsorbed fiux d 


2.7 x lO" 11 


1.8 x 10" 11 


1.2 x lO" 11 


1.4 x 10~ n 


Luminosity 


(3-l±i;g) x 10 33 D 2 5 








(0.5-2.5 keV) (erg/s) e 


(2.0ll; 6 ) x io 33 d 2 5 


(2.4+^ 2 ) x 10 34 D 2 8 


(2.7±i;l) x 10 34 D 2 8 


n e (cm -3 ) 


(0.11±0.03)(/D 3 . 5 )- 1/2 


(O.llt^tX/Ds.s)- 1 / 2 


(0.07i»;° 3 2 )(/D 4 . 8 )-V2 


(0.07±0.03)(/D 4 . 8 )- 1/2 


i f (kyr) 




(23i^)(//J 3 . 5 ) 1/2 


(30i 1 1 2 )(//3 4 . 8 ) 1/2 


(22l^)(/D 4 . 8 ) 1 /2 


Mass of X-ray 










emitting gas (Mq) 


(2.8 ± 1.5)/V2£)|/ 2 




{2.7t\l)f^D\^ 


(2.7±1.8)/ 1 / 2 D^ 8 2 



a (10 14 /4tt D 2 ) J n e nadV cm 5c " Elemental abundances arc relative to solar abundance. From X-ray extraction radius in Figures [2] and [5] 
.D3.5 and -D4.8 are distances in terms of 3.5 and 4.8 kpc, respectively. ' Flux in erg cm _2 s~ lc / is the volume filling factor/ Shock age 
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TABLE 2 

Catalogue of point sources surrounding SNR G85.4+0.7. Chandra positions are denoted by a CXO prefix 



# 


IAU Name 


OJ2000.0 


<5j2000.0 


Pos. Err. 
(arcsec) 


Counterpart 3 


Offset 
(arcsec) 


Optical 

A ray 
Flux Ratio 


1 


XMMU J205101. 6+452218 


20 51 01.593 


+45 22 18.47 


11.9 


1350-13030003 


1.3 


6.6 X 10 12 




CXO J205101. 4+452219 


20 51 01.392 


+45 22 18.60 


0.1 








2 


XMMU J205056. 8+452322 


20 50 56.750 


+45 23 22.18 


9.5 


1350-13027640 


2.2 


1.8 x 10 12 




CXO J205056.6+452320 


20 50 56.603 


+45 23 20.09 


0.6 








3 


XMMU J205043. 2+452213 


20 50 43.212 


+45 22 12.96 


11.2 


1350-13020622 


2.3 


6.1 x 10 12 




CXO J205043. 1+452215 


20 50 43.097 


+45 22 15.49 


0.8 








1 


XMMU J205034.8+451831 


20 50 34.811 


+45 18 31.47 


10.1 


1350-13016264 


2.9 


6.8 x 10 14 




CXO J205034. 9+451836 


20 50 34.945 


+45 18 35.83 


1.2 








5 


XMMU J205024. 5+452343 


20 50 24.523 


+45 23 42.95 


11.9 


1350-13011037 


2.7 


2.7 x 10 12 


6 


XMMU J205120. 9+451859 


20 51 20.925 


+45 18 58.73 


12.4 


J205120.79+451900.1(S) b 


2.4 


1.2 x 10 15 




CXO J205120. 8+451901 


20 51 20.802 


+45 19 01.11 


0.4 
















1350-13039935 


4.1 


2.6 x 10 


7 


XMMU J205127. 2+452025 


20 51 27.152 


+45 20 24.96 


11.6 


1350-13042959 


5.1 
1.7 


4.6 x 10 15 




/~t"V/~\ 7 on c 1 or 7 o i /i conic 

OaU JzU51z7.z+45zUz5 


on ci ot i co 

ZU 51 27. loo 


i a c on o c n 1 

+45 20 2b. (Jl 


0.6 








8 


XMMU J 205000. 7+452044 


on c n nn /^ct 

20 50 00.657 


i a c on a o n a 

+45 20 43.94 


12.2 


Toncnnn r* a i a con a /~ j n / o \ 

J 205000. 64+452046. 2(b) 


2.0 


1.4 X 10 












1 o c n innnoooi 

1350-12998831 


3.7 


"I A ^ , 1 nln 

1.4 X 10 


9 


"V TV /r i\ T T T ion ,1 ncn O I A COO /in 

aMMU JzU4y5y.o+45zo4y 


on a n en o nn 

20 49 59.800 


i a c oo a n nn 

+45 23 49.00 


14.5 


i o c n i onno \ n 1 

1350-12998401 


1.7 


a o w inl2 

4.o X ID 












1 o c o nonnonn/ T~> \ 

1353-0399326(B) 


2.1 


r* o , 1 nil 

6.8 X 10 


10 


nvA T on c 1 nn n i *i ro < /in 

CXO J205109. 0+452440 


on n nn noo 

20 51 09.022 


i /ic o a on oc 

+45 24 39.85 


1.4 


i of n i o n o o o o 

1350-13033336 


11.2 


o o ^ , 1 n1 .'^ 

2.3 X 10 ° 


11 


CXO J205056.6+451744 


on cn f i ' coo 

20 50 56.588 


I ^ r 1 7 A A 11 

+45 17 44.11 


1.0 


1 o c n 1 ono^T/inon 

1350-130274980 


6.3 


1.2 X 10 


12 


CXO J205HL9+452136 


20 51 11.876 


+45 21 36.26 


1.0 


1350-13035376 


1.7 


2.1 x 10 14 


13 


CXO J205036.9+451945 


20 50 36.913 


+45 19 45.00 


1.1 


1350-13017472 


7.8 


1.4 x 10 12 


14 


CXO J205138.4+452333 


20 51 38.370 


+45 23 33.43 


1.2 


1350-13049047 


6.0 


3.3 x 10 12 


15 


CXO J205H5.5+452207 


20 51 15.499 


+45 22 07.42 


1.8 


1353-0400904(B) 


6.2 


8.4 x 10 11 


16 


CXO J205103.6+451730 


20 51 03.641 


+45 17 30.00 


1.7 


1350-13031016 


3.8 


2.5 x 10 12 


17 


CXO J205149.3+452341 


20 51 49.264 


+45 23 40.59 


1.7 


1350-13054386 


2.4 


1.0 x 10 12 


18 


CXO J205140.0+452107 


20 51 40.002 


+45 21 06.64 


2.0 


1350-13049401 


5.6 


4.9 x 10 11 


19 


CXO J205143.3+452037 


20 51 43.338 


+45 20 37.04 


1.7 


1350-13051265 


3.2 


3.7 x 10 12 


20 


CXO J205135.2+451830 


20 51 35.164 


+45 18 30.27 


1.7 


1350-13047169 


1.5 


6.4 x 10 15 


21 


CXO J205126.4+451856 


20 51 26.420 


+45 18 56.06 


1.6 


1353-0401091 


7.2 


1.8 x 10 12 



a All designations are from the USNO-A2.0 catalogue JMonet et al.| I1998D unless otherwise noted. (S) denotes SKY2000 designation 
dMvers et al.||2002i ) c This row refers to the XMM-Newton position d This row refers to the Chandra position 
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TABLE 3 

Catalogue of point sources surrounding SNR G85.9— 0.6. Chandra positions are denoted by a CXO prefix 







IAU Name 


"J2000.0 


.12000.0 


r os. Err. 
(arcsec) 


Counterpart a 


(arcsec) 


Optical 

X — ray 
Flux Ratio 


1 


XMMU J205911.4+444729 


20 59 11.398 


+44 47 29.27 


11.3 


1347-0402736(B) b 


5.5 


1.1 X 10 19 


2 


XMMU J205950.5+445954 


20 59 50.504 


+44 59 53.58 


10.3 


HD 200102(I) C 


6.0 


3.7 x 10 15 


3 


XMMU J205944.3+450741 


20 59 44.261 


+45 07 40.75 


11.7 


1350-13269810 


3.7 


1.1 x 10 14 


1 


XMMU J205857.0+450348 


20 58 56.993 


+45 03 48.25 


11.0 


1350- 0403688(B) 


3.7 


1.0 x 10 13 


5 


XMMU J210003. 3+450342 


21 00 03.280 


+45 03 41.76 


10.4 


1350-13276811 


1.3 


2.5 x 10 12 


6 


XMMU J210023. 1+445435 


21 00 23.149 


+44 54 34.98 


11.1 


21002305+4454364(M) d 


1.4 




7 


XMMU J210000. 1+445631 


r> i r\r\ nn ntT 

21 00 00.057 


i a a tic on n a 

+44 56 30.94 


9.3 


21000013+445o3o4(M) 


4.9 




8 


XMMU J205917.0+445809 


, 1 / 1 ( 1 1 r\r\f\ 

20 59 17.009 


1 A A P"0 no *~7 ~1 

+44 58 08.71 


10.5 


1349- 0404180(B) 


1.4 


4.2 X 10 




CXO J205917.3+445821 


on Kfl 1 *7 one 

20 59 1 / .308 


i/i/i £ q on on 
+44 58 20.90 


0.5 


12.2 




e 
f 












134y-U4U4iyi(r> ) 


11.7 


2.7 X 10 14 


f 












1275-14437760 


1.1 

12.3 


5.5 x 10" 


9 


CXO J205915. 7+445858 


20 59 15.715 


+44 58 57.64 


1.1 


1275-14437093 


0.7 
9.6 


6.2 x 10 14 


10 


CXO J205849. 8+445725 


20 58 49.768 


+44 57 25.00 


0.8 


1275-14424519 


2.7 


8.2 x 10 13 


11 


CXO J205839.9+445132 


20 58 39.924 


+44 51 31.67 


1.3 


1275-14419711 


5.1 


1.8 x 10 15 


12 


CXO J205835. 3+444640 


20 58 35.349 


+44 46 40.21 


1.2 


1347-0402278(B) 


1.1 


8.4 x 10 11 


13 


CXO J205825. 4+444940 


20 58 25.398 


+44 49 39.63 


1.7 


1348-0403023(B) 


7.1 


6.1 x 10 14 


11 


CXO J205823. 9+444621 


20 58 23.867 


+44 46 21.35 


1.7 


1347-0402123(B) 


3.2 


1.8 x 10 12 














20582367+4446195(M) 


1.8 


2.9 x 10 13 


15 


CXO 


J205930.8+445729 


20 59 30.760 


+44 57 29.19 


1.1 


1349-0404462(B) 


4.0 


16 


CXO 


J205922. 2+445416 


20 59 22.191 


+44 54 15.95 


1.5 


1349-0404292(B) 


3.1 


1.8 x 10 13 














20592216+4454176(M) 


1.7 


2.6 x 10 12 


17 


CXO 


J205921. 9+445818 


20 59 21.928 


+44 58 18.35 


1.1 


1275-14439638 


8.6 














20592206+4458188(M) 


1.6 


7.0 x 10 13 


18 


CXO 


J205920.0+445321 


20 59 19.993 


+44 53 21.12 


1.3 


1275-14439333 


7.0 














20592022+4453178(M) 


4.1 




19 


CXO 


J205912. 6+445033 


20 59 12.619 


+44 50 33.10 


1.1 


20591251+4450356(M) 




3.5 x 10 13 


20 


CXO 


J205904. 1+445952 


20 59 04.133 


+44 59 52.38 


1.1 


1275-14431663 


7.8 


21 


CXO 


J205857. 7+450008 


20 58 57.717 


+45 00 07.51 


1.2 








22 


CXO 


J205851. 5+450050 


20 58 51.450 


+45 00 50.04 


1.3 


1350-0403625(B) 


3.5 


2.4 x 10 12 


23 


CXO 


J205848. 7+445920 


20 58 48.681 


+44 59 19.84 


1.2 


1275-14423599 


8.4 


2.3 x 10 13 


24 


CXO 


J205841. 1+445260 


20 58 41.091 


+44 52 59.50 


1.1 








25 


CXO 


J205837.4+445510 


20 58 37.389 


+44 55 09.86 


1.1 


1275-14418321 


9.1 


1.0 x 10 13 














20583749+4455092(M) 


1.3 


3.2 x 10 12 


26 


CXO 


J205834.6+445330 


20 58 34.649 


+44 53 30.33 


1.4 


1348-0403158(B) 




27 


CXO 


J205814.7+444606 


20 58 14.652 


+44 46 05.70 


2.1 









a All designations are from the USNO-A2.0 catalogue (Monet ct al. 1998) unless otherwise noted. (B) denotes USNO-B1.0 designation jMonet et al.| 
I2003l) c (I) denotes designation in lGuarinosl <|1992D d (M) denotes 2MASS designation JSkrutskie et al.|[2006h e This row refers to the XMM-Newton 
position^ This row refers to the Chandra position 
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TABLE 4 

Parameters for fits of XMM-Newton PN, MOS and Chandra 3 - point source spectra to an absorbed power law. 



# N u (10 22 cm~ 2 ) T PL Normalization 13 \l 0) 0.5-2.0 keV 2.0-10.0 keV Hardness ratio c X-ray 

Counts Counts S/N 

G85.4+0.7 



1 


50+ - 25 
u - ou -0.18 


z - 4 '-0.47 


4.75t?J 


1.20 (37) 


210 


56 


0.27 


3.8 


2 


94+ - 27 




11.2l» 2 8 


0.93 (14) 


122 
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a Chandra spectra are included in fit for sources 1, 2,3, 4, 6, and 7 of G85.4+0.7 and source 8 of G85.9— 0.6.^ 10 5 photons keV 1 cm 2 s 1 
at 1 kcV c Hardness ratio is the ratio of hard counts to total (hard+soft) counts 
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Fig. 1. — 0.5—2.5 keV XMM-Newton X-ray mosaic image of G85.4+0.7 with the highest few radio contours overlaid in blue. The X-ray 
image is Gaussian smoothed with a radius of 3 pixels. 
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Fig. 2. — PN and MOS mosaic images of G85.4+0.7. Top panel: 0, 5-2. 5 keV image, showing point sources (1-9) and the diffuse emission 
extraction region, as well as the region used for the central blob (see i|7.3l l. The background region for the diffuse spectra is also shown as a 
dashed circle. The image is smoothed as in Figure[l] Middle panel: 2.5-10.0 keV image, similarly smoothed and with the same contrast as 
in the top panel. Bottom panel: 0.5—2.5 keV image, Gaussian smoothed with a radius of 3 pixels and logarithmically scaled and adjusted 
in contrast to emphasi ze the diffuse emissi on and overlaid with contours of the 0.5-2.5 keV emission smoothed similarly to the ROSAT 
contours in Figure 6 of lKothes etaLl [EfflOlT ). 
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Fig. 3. — Chandra image of G85.4+0.7 from the S2 and S3 chips, Gaussian smoothed with a radius of 5 pixels, showing (in black) point 
source extraction regions and corresponding regions (in blue) for XMM-Newton data for comparison. The large circle indicates the diffuse 
region in the top panel of Figure[2]and the background regions are shown as dashed circles. (See the online version for colour figure.) 
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Fig. 4. — 0.5—2.5 keV XMM-Newton X-ray mosaic image of G85.9— 0.6 with the highest few radio contours overlaid in blue. The X-ray 
image is Gaussian smoothed with a radius of 3 pixels. 
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Fig. 5. — PN and MOS mosaic images of G85.9— 0.6. Top panel: 0. 5—2. 5 keV image, showing point sources (1-8) and the diffuse emission 
extraction region, as well as the region used for the central blob (see i|7.4p . The background region for the diffuse spectra is also shown as a 
dashed circle. The image is smoothed as in Figure[4] Middle panel: 2.5-10.0 keV image, similarly smoothed and with the same contrast as 
in the top panel. Bottom panel: 0.5—2.5 keV image, Gaussian smoothed with a radius of 3 pixels and logarithmically scaled and adjusted 
in contrast to emphasi ze the diffuse emissi on and overlaid with contours of the 0.5-2.5 keV emission smoothed similarly to the ROSAT 
contours in Figure 6 of Kothcs ct al.l 1)200 lh . 




Fig. 6. — Chandra image of G85.9— 0.6 from the S2 and S3 chips, Gaussian smoothed with a radius of 5 pixels. Extraction regions for 
sources 1 and 8 from the XMM-Newton image are shown as blue circles. (See the online version for colour figure.) 
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Fig. 7. — XMM-Newton PN background subtracted 0.5—8.0 keV spectra of G85.4+0.7, using backgrounds from the observation (black) 
and from blank sky event files (grey). The oversubtraction by the blank sky background is evident. The spectrum using the observation 
background is clearly background dominated above ~2.5 keV. 
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Fig. 8.— XMM-Newton PN and MOS 0.5-2.5 keV spectra of G85. 4+0.7, fit to an absorbed VPSHOCK. The background spectrum for 
PN was extracted in the same way as for the sp ectru m shown in Figure |91 and the background for the MOS instruments are extracted from 
the XMM-Newton ESAS package described in i|4.1l The fitted parameters are given in the third column of Table [T] The PN is shown in 
black, and the MOS 1 and 2 are shown in blue and red res pect ively. The bottom panel displays the residuals in units of a. Instrumental 
lines in the MOS spectra appear at 1.49 and 1.75 keV (see ^14. II for details). The line above 1.75 keV is an elemental line which appears in 
spectra from all instruments. 
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Fig. 9. — XMM-Newton PN and MOS spectra of G85.4+0.7, fit to an absorbed VPSHOCK. The background spectra used here were 
extracted from the region shown in Figure [2] The fitted parameters are given in the second column of Table [T] The PN is shown in black, 
and the MOS 1 and 2 are shown in blue and red respectively. The bottom panel displays the residuals in units of a. An instrumental line 
in the MOS spectra appears at 1.75 keV (see ^14. II for details). 
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Fig. 10.— XMM-Newton PN and MOS 0.5-2.5 keV spectra of G85. 9-0.6, fit to an absorbed VPSHOCK. The background spectrum for 
PN was extracted in the same way as for the spect rum shown in Figure 1111 and the background for the MOS instruments are extracted 
from the XMM-Newton ESAS package described in EI4.1I The fitted parameters are given in the third column of Table [T] The PN is shown 
in black, and the MOS 1 and 2 are shown in blue and red re spec tively. The bottom panel displays the residuals in units of a. Instrumental 
lines in the MOS spectra appear at 1.49 and 1.75 keV (see ^14. II for details). The line above 1.75 keV is an elemental line which appears in 
spectra from all instruments. 



XMM-Newton, Chandra, & CGPS observations of G85.4+0.7 and G85.9-0.6 



25 




Fig. 11.— XMM-Newton PN and MOS 0.5-2.5 keV spectra of G85.9-0.6, fit to an absorbed VPSHOCK. The background spectra 
used here were extracted from the region shown in Figure [5] The fitted parameters are given in the fourth column of Table [T] The PN is 
shown in black, and the MOS 1 and 2 are shown in blue and red respectively. The bottom panel displays the residuals in units of a. An 
instrumental line in the MOS spectra appears at 1.75 keV (see i|4,ll for details). The excess counts in the MOS spectra around 1.4 and 0.85 
keV are most likely due to instrumental effects. 
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Fig. 12. — XMM-Newton PN, MOS, and Chandra 0.5-5 keV spectra of G85.4+0.7 source 1, fit to an absorbed power law. The fitted 
parameters are given in Table [4] The PN is shown in solid black, the MOS 1 and 2 are shown in dashed blue and solid red respectively, 
and the Chandra spectrum is shown as dotted black. The bottom panel displays the residuals in units of a. 
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Fig. 13. — Integrated foreground hydrogen column density profiles as a function of radial velocity for the two SNRs G85.4+0.7 (left: solid 
and dotted line) and G85.9— 0.6 (right: solid line). The dashed lines represent the upper and lower limits for the absorbing H I column 
density determined from the X-ray spectra. The dotted-dashed line in the diagram for G85.4+0.7 represents the radial velocity of the 
related stellar wind bubble determined by Kothes et al. (2001). 
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Fig . 14.— The velocity field of the H I gas towards £ =86.0°, 6 = -0.6°, from the H I modelling method of IFoster fc Mac Williams] 
(2006). SNRs G85.4+0.7 & G85.9— 0.7 are marked with triangles showing their respective distances of 3.5 and 4.8 kpc. The Perseus Spiral 
arm's peak density is at 5.0 kpc and the spiral shock that precedes the arm is at 4. 8 kpc. The solid dark line is the underlying circular 
velocity field fitted in this direction, and is a power-law in galactoccntric distance (sec Foster & MacWilliams 2006) with a gradient similar 
to a Schmidt rotation curve. 
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Fig. 15. — 0.5—2.0 keV XMM-Newton X-ray mosaic image of G85.4+0.7 with radio contours overlaid in blue. The X-ray image has had 
the point sources removed and is Gaussian smoothed with a radius of 1' to match the radio resolution, in order to compare the radio and 
diffuse X-ray emission regions. The radio contours were chosen to approximate the appearance of the top image in Figure 2 in lKothes et al.l 
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Fig. 16. — 0.5—2.5 keV XMM-Newton X-ray mosaic image of G85.9— 0.6 with radio contours overlaid in blue. The X-ray image has had 
the point sources removed and is Gaussian smoothed with a radius of 1' to match the radio resolution, in order to compare the radio and 
diffuse X-ray emission regions. The radio contours were chosen to approximate the appearance of the top image in Figure 4 in lKothes et al.l 



